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ABSTRACT 


The  coupling  of  magnetohydrodynamic  waves  to  electromagnetic 
and  sound  waves  at  a  plane  interface  between  a  plasma  and  a  neutral 
gas  is  considered  for  arbitrary  direction  of  the  incident  wave  vec¬ 
tor  and  for  the  constant  external  magnetic  field  lying  either  in  the 
plane  of  incidence  or  perpendicular  to  it.  The  relations  between  the 
angle  of  incidence,  reflection,  and  transmission,  the  ratio  of  the 
field  amplitude  transmitted  or  reflected  to  the  corresponding  ampli¬ 
tude  incident,  and  the  energy  coupling  coefficients  are  calculated. 

An  incident  fast  magnetoacoustic  mode  generates  a  sound  wave 
whose  wave  propagation  vector  makes  an  angle  with  the  normal  to  the 
boundary  of  the  order  of  the  sound  speed  to  Alfven  speed,  when  this 
ratio  is  small.  In  this  case,  the  energy  coupling  coefficient  of 
sound  to  the  fast  magnetoacoustic  mode  is  also  of  the  order  of  this 
ratio.  An  incident  slow  magnetoacoustic  mode  generates  a  sound  wave 
with  energy  coupling  coefficient  of  the  order  of  unity,  and  an  elec¬ 
tromagnetic  wave  with  energy  coupling  coefficient  of  the  order  of 
the  ratio,  speed  of  sound  to  speed  of  light. 
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1 .  INTRODUCTION 


Although  magnetohydrodynamic  wave  propagation  has  been  thoroughly 
treated  for  the  case  of  a  homogeneous  infinite  medium,^  in  many  cases  of  prac¬ 
tical  significance  the  results  are  not  applicable  since  the  medium  is  inhomo¬ 
geneous  due  to  spatial  changes  in  the  ionization  density.  The  present  paper 
i s  concerned  with  wave  propagation  in  a  medium  with  a  spatial  discontinuity 
in  the  ionization  density.  In  particular,  this  discontinuity  will  be  taken 
so  that  the  gaseous  medium  changes  from  one  which  is  fully  ionized  to  one 
which  1 s  totally  un-ionized.  The  ionized  medium  can  support  three  distinct 
magnetohydrodynamic  modes  (the  Alfven,  the  fast  magnetoacoustic,  and  the  slow 
magnetoacoustic),  whereas  the  un-ionized  medium  can  support  electromagnetic 
waves  and  sound  waves.  The  coupling  of  the  magnetohydrodynamic  modes  to  the 
electromagnetic  and  acoustic  modes  at  a  plane  plasma-neutral  gas  interface 
will  be  discussed  in  the  present  paper. 

This  problem  is  similar  to  two  others  which  also  involve  discontin¬ 
uities  in  the  gaseous  medium.  All  three  may  be  described  in  terms  of  the  ion¬ 
ization  concentration  on  each  side  of  the  interface: 

A.  the  media  on  both  sides  of  the  interface  are  fully  ionized  (plasma-plasma 
problem) ; 

B.  the  medium  from  which  the  waves  are  approaching  the  interface  is  totally 
un-ionized,  whereas  the  other  medium  is  fully  ionized  (neutral  gas-plasma 
problem) ; 

C.  the  medium  from  which  the  waves  are  approaching  the  interface  is  fully 
ionized,  whereas  the  other  medium  is  totally  un-ionized  (plasma-neutral  gas 
problem) . 
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A.  Plasma-Plasma  Problem 

Any  of  the  three  magnetohydrodynamic  modes  (Alfven,  fast  magneto¬ 
acoustic,  and  slow  magnetoacoustic)  may  be  incident,  and  in  the  most  general 

2 

case  all  three  modes  can  be  reflected  or  transmitted.  Ferraro  considers 
reflection  and  transmission  of  Alfven  waves  in  incompressible  media  for  the 
magnetic  field  of  the  incident  wave  parallel  to  the  interface.  Roberts^ 
extends  Ferraro's  work  to  arbitrary  orientations  of  the  magnetic  field  of 
the  incident  wave.  Simon\  Williams^,  and  Raju  and  Verma^  consider  reflec¬ 
tion  and  transmission  of  all  three  modes  in  compressible  media  and  examine 
in  detail  the  case  where  the  external  magnetic  field  is  normal  to  the  in- 

7 

terface.  Talwar  extends  this  analysis  by  assuming  the  sound  velocity  is 

/  O 

much  smaller  than  the  Alfven  velocity.  Fejer°  examines  dispersion  relations 
and  Snell's  laws  when  the  two  plasma  regions  have  a  relative  velocity  paral¬ 
lel  to  the  interface.  He  chooses  the  external  magnetic  field  parallel  to 
the  interface.  Frieman  and  Kulsrud^  calculate  amplitude  relations  for  all 
three  modes  when  the  external  magnetic  field  is  parallel  to  the  interface 
and  lies  in  the  incident  plane.  Pr idmore-Brown^  gives  a  graphical  presen¬ 
tation  of  Snell's  laws  to  determine  the  nature  and  disposition  of  reflected 
and  refracted  waves.  Bazer^  calculates  amplitude  relations  in  detail  for 
the  external  magnetic  field  in  the  incident  plane. 

B.  Neutral  Gas-Plasma  Problem 

The  two  electromagnetic  modes  and  the  sound  mode  may  be  incident 

on  the  interface  and  reflected  from  it,  while  all  three  magnetohydrodynamic 

12 

modes  can  be  transmitted.  Turcotte  and  Schubert  calculate  the  coupling 
of  an  Incident  electromagnetic  wave  with  magnetohydrodynamic  waves.  They 
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consider  the  cases  where  the  external  magnetic  field  is  parallel  to  the  In¬ 
terface,  the  magnetic  field  of  the  incident  wave  is  parallel  to  the  external 
magnetic  field,  and  the  incident  wave  vector  is  either  normal  to  or  parallel 
to  the  interface.  Kornhauser^  considers  in  greater  detail  the  case  where 
the  external  magnetic  field  is  parallel  to  the  interface  and  the  incident 
wave  vector  is  normal  to  the  interface.  Kontorovich  and  Glutsyuk^  calcu¬ 
late  the  coupling  of  an  incident  sound  wave  with  magnetohydro dynamic  waves 
/ 

for  the  Alfven  velocity  much  less  than  the  sound  velocity. 

C.  Plasma-Neutral  Gas  Problem 

Any  of  the  three  magnetohydrodynamic  modes  can  be  incident  on  and 
reflected  from  the  interface,  while  acoustic  waves  and  electromagnetic  waves 
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of  both  polarizations  can  be  transmitted.  Kahalas  considers  the  case  in 

which  the  sound  velocity  is  zero  and  both  the  incident  wave  vector  and  the 

external  magnetic  field  are  normal  to  the  interface,  Ullah  and  Kahalas^, 
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and  Kahalas  '  discuss  the  case  in  which  the  sound  velocity  is  zero,  the  in¬ 
cident  wave  vector  has  arbitrary  orientation,  and  the  external  magnetic  field 
lies  either  in  the  plane  of  incidence  or  perpendicular  to  it.  In  the  above 
analyses  only  the  electromagnetic  waves  are  transmitted  since  with  sound  vel¬ 
ocity  zero  sound  waves  do  not  propagate. 

This  paper  is  an  extension  of  the  latter  work  to  the  case  of  non¬ 
zero  sound  velocity.  The  incident  wave  vector  will  be  allowed  to  have  an 
arbitrary  orientation  with  respect  to  the  interface,  but  the  constant  ex¬ 
ternal  magnetic  field  will  be  constrained  to  lie  either  in  the  plane  of  in¬ 
cidence  or  perpendicular  to  it. 


The  mathematical  model  analyzed  in  this  paper  has  features  in  com- 


mon  with  the  physical  situation  which  describes  the  propagation  of  geomag¬ 
netic  disturbances  down  through  and  below  the  ionosphere.  The  ionosphere 
and  atmosphere  are  considered  here  as  two  distinct  regions,  with  the  ionos¬ 
phere  as  a  homogeneous  plasma,  the  atmosphere  as  a  homogeneous  neutral  gas, 
and  the  transition  between  the  two  as  a  plane  surface. 

Within  the  ionosphere  and  the  atmosphere,  the  sound  velocity  is 
much  less  than  the  Alfven  velocity;  consequently,  this  fact  will  be  used 
to  simplify  some  of  the  equations. 

The  physical  model  discussed  here  is  lacking  because  the  external 
magnetic  field  is  assumed  to  be  uniform  whereas  the  earth's  field  is  not, 
collisions  between  positive  ions  and  neutral  molecules  in  the  lower  ionos¬ 
phere  have  not  been  taken  into  account,  the  ionospheric  boundary  is  not 
necessarily  as  abrupt  as  assumed,  and  the  effect  of  the  conducting  earth 

itself  has  not  been  included.  Numerical  studies  have  been  made  by  Francis, 

1 8 

Karplus,  and  Dragt  on  the  propagation  of  magnetohydrodynamlc  waves  through 
the  lower  ionosphere  incorporating  the  positive  ion-neutral  particle  colli¬ 
sions,  the  effect  of  the  earth,  and  the  gradual  change  of  ionospheric  param¬ 
eters.  The  present  analytical  approach  should  be  considered  as  a  supplement 
to  their  work. 

In  Section  2,  a  derivation  is  given  of  the  boundary  conditions  re¬ 
lating  field  quantities  across  the  fluid  interface.  These  boundary  condi¬ 
tions  are  valid  for  the  plasma-plasma,  neutral  gas -plasma,  and  plasma-neu¬ 
tral  gas  problems.  Section  3  gives  an  overall  statement  of  the  problem. 
Section  4  discusses  the  properties  of  magnetohydrodynamic  waves.  In  par¬ 
ticular,  the  dispersion  relations,  relation  between  field  quantities,  and 
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energy  transport  vector  are  exhibited.  Section  5  discusses  the  same  topics 
for  the  electromagnetic  and  acoustic  waves  propagating  in  the  neutral  gas. 
Section  6  treats  the  relations  between  the  angles  of  incidence,  reflection, 
and  transmission  and  the  relations  between  the  field  quantities  on  both 
sides  of  the  interface  for  the  case  in  which  the  constant  external  magnetic 
field  lies  normal  to  the  plane  of  incidence.  Section  7  similarly  discusses 
the  case  in  which  the  constant  external  magnetic  field  lies  in  the  plane 
of  incidence. 

2.  BOUNDARY  CONDITIONS 

The  boundary  conditions  relating  the  field  quantities  on  two  sides 

of  a  discontinuity  have  been  derived  by  several  authors  using  somewhat  dif- 

19,20 

ferent  methods.  '  ’  The  present  derivation  involves  an  integration  of  the 
conservation  equations  across  the  boundary.  The  boundary  conditions  which 
result  are  equally  valid  for  the  plasma-plasma,  plasma-neutral  gas,  and  neu¬ 
tral  gas-plasma  cases. 

In  order  to  derive  the  boundary  conditions  at  a  plasma-neutral  gas 

interface,  it  is  important  to  ensure  that  the  equations  used  are  valid  in 

both  the  plasma  and  the  neutral  gas.  This  means,  for  example,  that  the3E/3t 

term  (needed  to  describe  electromagnetic  waves)  must  be  included  in  the  ?xB^ 

Maxwell  equation,  although  it  is  often  left  out  of  the  magnetohydrodynamic 

equations.  Similarly,  Ohm's  law  is  written  with  a  finite  conductivity, 

since  the  variation  of  the  conductivity  (from  infinity  in  the  plasma  to  zero 

21 

in  the  neutral  gas)  must  be  included.  The  equations  used  are 

pfav/et)  -  -p(v-v)^  -  vp  +  (  JxB)  /  c  +  pi  (2.1) 
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(2.2) 


and 


(bp/bt )  +  v- Cpv)  =  0 

vxE  =  -(l/c)(dB/dt)  ,  (2.3) 

vxB  =  (l/c)(8E/6t)  +  Uflj/c  ,  (2.M 

M. 

vi-0  ,  (2.5) 

v-E  =  Uitp  ,  (2.6) 

•MJ  c 

J  =  <r(E  +  VXB/C)  ,  (2.T) 

vp  =  vp  ,  (2.8) 


where  E,  B,  v,  p,  j,  p,  p  c,  s,  and  a  are  total  electric  field,  magnetic 

^  Wv  w  c 

field,  fluid  velocity,  pressure,  electric  current  density,  mass  density,  elec 
trie  charge  density,  velocity  of  light,  velocity  of  sound,  and  electrical  con 
ductivity,  respectively.  A  bar  over  a  field  quantity  denotes  a  total  unlin¬ 
earized  quantity,  a  subscript  zero  denotes  constant  quantities,  and  field 
quantities  without  bars  are  used  to  denote  quantities  of  small  amplitude. 

These  equations  may  be  rearranged  to  express  the  conservation  of  mo- 

19 

mentum 


(d/dt)  +  (l/Unc)  Ex_B 


Z  «/6x. 

J  J  J 


(2.9) 


and  conservation  of  energy 

(b/bt)  [cpv2/2)  +  (r  +  ‘r)/3jTj  +  v-T  -  p(v.v)  +  j-/a  -  (v-I)(vE)Ajt  -  0  (2.10) 
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where 
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system.  Then  Eqs.  (2.13)  —  (2.16)  may  be  written  to  first  order  in  the  per¬ 

turbation  as 


A  A 

nxE.  =  nxE0  , 

■'*“  •«*.  vw 

(2.17) 

A  A 

(2.18) 

v  *n  =  v2'n 

A  A 

Bi"n  =  B2‘n  • 

(2.19) 

and 


[P1  +  (B0-B1)/4*j2  -(l/4it)|j^(B0-£)  +  W£)] 

=  ^p2  +  (Bo-32)/4nJ^  -  (l/^Tt)  j^>( BQ'n)  +  B0(B2.^  .  (2.20) 

A  A  A  A  .  .  .  . 

With  the  vector  identity  B  =  (B*n)n  -  nx(nxB)  along  with  Eq .  (2.19)»  Eq.  (2. 20} 

<*w  «»>  -«w  ^  -r>-- 

may  be  separated  into  its  normal  and  tangential  parts,  giving 


P1  +  (Bo'Bl)A,t  =  p2  +  (B0-B2)/4n 


and 


A  A 


A  A 

(B  -n)  nxB.  =  (B  -n)  nxB~ 
o  1  o'..  2 


(2.21) 

(2.22) 


» 

There  exist  two  separate  cases  depending  on  whether  BQ*n  »  0  or  BQ*n  r  0* 

When  BQ.n  -  0,  the  boundary  conditions  are  Eqs.  (2.17),  (2.l8),  and 

«***» 

(2.21).  A  discontinuity  in  B,  resulting  from  a  surface  current  density,  J*, 


is  determined  by  nx(B2  -  B^)  =  4jtj*/c. 


When  B0-n  =  0,  the  boundary  conditions  reduce  to  Eqs.  (2.17)  and 


(2. l8)  along  with 


and 


(2.23) 

(2.24) 
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It  will  be  shown  now  that  Eq .  (2.10)  is  consistent  with  the  boundary 
conditions  to  first  and  second  order  in  the  perturbation.  To  first  order  Eq. 
(2.10)  is  the  scalar  product  of  Eq.  (2.3)  with  Bq.  The  resulting  boundary  con 
dition  duplicates  the  scalar  product  of  Eq.  (2.17)  with  BQ.  To  second  order, 
Eq.  (2.10)  becomes 


■aw/at  +  v  •  T  +  j  2/<r  =  0  , 

(2.2?) 

where 

T  =  pv  +  c(ExB)Ajt  , 

(2.26) 

and 

W  2  pQv2/2  +  (E2  +  B2)/8n  +  s2p2/2pQ 

The  integration  of  Eq.  (2.25)  over  a  Gaussian  pillbox  of  base  area  &A  and  vol¬ 
ume  Sm  extending  through  the  surface  gives 


A  A 

T  *n  *  Tg-n  +  A 


(2.27) 


where  A  3  him  (l/SA).Cl /C  dV(j2/4r").  For  the  case  B  -n  ^  0  there  is  no  surface 
5A-»0  u  - 

current  density  and  A  =  0.  Then  Eq.  (2.27)  is  seen  to  be  consistent  with  Eqs, 
(2.17),  (2.18),  (2.23),  and  (2.2b).  For  the  case  B0'n  =  0,  there  is  a  sur- 
face  current  density.  A  is  calculated  directly  from  j2/o"  ■  j-(E  +  vxBQ/c) 
evaluated  over  the  surface.  The  surface  current  density  is  by  definition  tan¬ 
gential  to  the  interface  and  (E  +  vxB0/c)  .  .  is  continuous  across  the 

-1^-  tangential 

boundary  so  that  A  is  given  by  A  ■  j*-(E^+  v^xBQ/c).  The  term  in  parentheses 

**  -w. 

vanishes  when  evaluated  in  the  plasma  medium  where  the  conductivity  is  infln- 

A  A  /s 

ite  so  that  A  “  0  and  Eq.  (2.27)  becomes  T, -n  ■  Tp-n.  If  (Tp  -  T, )n  is  cal- 
culated  directly  from  the  definition,  Eq,  (2.26),  using  Eqs.  (2,17)i  (2.18) 
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and  (2.21),  one  gets 

(T2  -  Tx)-n  =(c/4n)(B2  -  B1)  •fnxE1  -  (n-v^cjB  "1  . 

Again,  the  bracketed  term  vanishes  when  the  conductivity  is  infinite  so  that 
Eq.  (2,10)  is  consistent  with  the  boundary  conditions  to  first  and  second 
order  in  the  perturbation. 

3.  DESCRIPTION  OF  PROBLEM 

In  the  ionized  medium,  the  equations  governing  extremely  low  fre¬ 
quency  wave  motion  (i.e.,  frequencies  less  than  the  ion  cyclotron  frequency) 
are  given  by  Eqs.  (2.1)  -  (2.8)  written  to  first  order  in  the  perturbation 
and  are 


p0(av/*t)  =  -vp  +  JxBp/c 

(3.D 

a  p/at  +  pQ**v  =  0  , 

(3.2) 

vxE  =  - ( l/c)3B/8t  , 

(3.3) 

vxB  =  4jrj/c  , 

(3.4) 

O 

II 

«| 

> 

(3.5) 

▼  •E  =  4jtp6  , 

(3.6) 

E  =  -vxB  /c  , 

—  -  *S> 

(3.7) 

„  2 

7p  ■  8  7p  , 

(3.8) 

where  infinite  conductivity  has  been  assumed  and  the*E/®t  term  in  Eq .  (2.4) 
has  been  neglected. 
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In  the  neutral  gas  the  conductivity  is  assumed  to  be  zero,  so  that 
the  electric  charge  and  current  densities  are  zero.  The  electromagnetic  equa¬ 
tions  are  then  uncoupled  from  the  fluid  flow  equations,  and  the  system  is  de¬ 
scribed  by 


Pofcv/aO  = 

? 

-sc  vp 

(3.9) 

ap/»t  +  p0v.v^  = 

0 

(3.10) 

vxE  = 

-( 1/c)  0B/dt  , 

(3.11) 

vxB  = 

(i/c)0E/at  , 

(3.12) 

V-  B  = 

0 

(3.13) 

v-E  = 

0 

(3.1*+) 

Since  Eqs.  (3.1)  -  (3.1*0  are  linear  and  homogeneous,  the  solutions  may  be  ob¬ 
tained  by  assuming  that  the  field  quantities  vary  as  exp  [i(k.r  -  «*)]  where  to 
is  the  angular  frequency  andjk  is  the  wave  vector.  Then  Eqs.  (3.1)  *  (3.8) 
describe  magnetohydrodynamic  waves  of  three  modes,  Eqs.  (3.9)  and  (3.10)  de¬ 
scribe  a  sound  wave,  and  Eqs.  (3.11)  -  (3.1*+)  describe  electromagnetic  waves 
with  two  independent  polarizations. 

The  geometry  is  chosen  so  that  the  x-y  plane  at  z  ■=  0  is  the  plane 
of  the  interface,  and  the  normal,  *n,  lies  along  the  positive  z  axis  and  points 
from  the  plasma  into  the  neutral  gas.  Without  loss  of  generality  k  will  be 
taken  to  lie  in  the  y-z  plane  which  is  the  plane  of  incidence.  Let  6  be  the 
angle  made  by  k  with  the  z  axis.  Then  k  may  be  written 


k  =  k(0,  sin  9,  cos  9) 


(3.15) 


The  calculations  will  be  made  for  two  cases:  BQ  perpendicular  to 

AIM 

the  plane  of  incidence  and  Bq  in  this  plane.  In  the  first  case,  only  the  fast 
magnetoacoustic  mode  propagates,  and  in  the  second  case  the  Alfven  mode  is  un¬ 
coupled  from  the  magnetoacoustic  modes. 

The  geometry  is  shown  in  Fig.  1  for  the  case  in  which  BQ  lies  in  the 
plane  of  incidence.  Superscripts  i,  r,  and  t  refer  to  the  incident,  reflected, 
and  transmitted  waves,  respectively.  There  is,  in  general,  more  than  one  mode 
reflected  and  transmitted  (  for  a  given  incident  mode)  but  this  is  not  shown 
in  the  figure.  The  ratios  of  the  transmitted  to  incident  fields  and  the  ra¬ 
tios  of  energy  transmitted  in  any  mode  to  the  incident  energy  will  be  cal¬ 
culated.  The  energy  coupling  coefficient  is  defined  to  be  the  latter  ratio, 

t  A  f  A 

<  T  -n  >/<  T  -n  >,  where  the  brackets  denote  a  time  average. 

>*V  **  AMI 

It  may  happen  that  for  certain  directions  of  k1,  some  kr  and  kfc  will 

<MA>  AdA" 

have  imaginary  z  components.  This  corresponds  to  the  physical  fact  that  these 
reflected  and  transmitted  waves  are  exponentially  attenuated  in  the  z  direc¬ 
tion.  The  attenuation  is  of  the  form  exp(-|z|/d)  where 

d=l/|kz|  .  (5.16) 

/  2  yk 

Now  k  *  k( 1  -  sin  6)  .  Whenever  sin  6  »  1,  then 
z 

d*  (k  (sir.  ©|)_1  =  |ky  I"1  .  (3. IT) 
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4.  PROPERTIES  OF  MAGNETOHYDRODYNAMIC  WAVES 


Equations  (3-1)  *  (3.8)  give 


a2v/at2  = 


s2v(  V 


v) 


(4.1) 


where  the  caret  denotes  a  unit  vector,  s  is  the  sound  velocity,  and  a  is  the 

/  L 

Alfven  velocity,  B0/(4itpo)  .  Equation  (4.1)  can  be  written 


£  M  v  =  0  ,  (4.2) 

j 

where 

Mjj  a  |^2  -  ^(Bq.Jc)2]^  -  (s2  +  a2)k1kJ  +  a2(Eo-^k)  |jc1Bo^  +  B0ikj  1*(4.5) 

The  condition  for  a  non-trivial  solution  of  Eq.  (4,2)  is  that  det  M  vanish; 
this  gives  the  dispersion  relations  for  three  modes 


U) 


» a|y  Bg 


(4.4) 


where  subscripts  is  used  to  denote  the  Alfven  mode, 


and 


=  k+a 

1  +  (s2/2a2)(k  xB  )2  +  0(s/a)M 

(4.5) 

used  i 

to  denote  the  fast  wave  of  the  magnetoacoustic 

mode, 

k-„ 

|  jl  -  (s2/2a2)(k_xB0)2  +  0(s/a)^j 

(4.6) 

where  subscript  -  is  used  to  denote  the  slow  wave  of  the  magnetoacoustic  mode. 

2 

The  dispersion  relations  are  exhibited  as  a  power  series  in  (s/a)  for  appli¬ 
cation  to  the  geophysical  problem  where  s/a  «  1. 
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The  eigenvectors  associated  with  each  mode  are  determined  from  Eq. 
(4.2)  with  the  aid  of  the  dispersion  relations.  The  results  are  that  is 
perpendicular  to  both  and  Bn,  and  v+  and  v_  lie  in  the  plane  containing 

w 

B0  and  the  corresponding  wave  vector. 

When  Bq  is  normal  to  the  plane  of  incidence,  i.e.,  BQ  »  BQ(l,  0>  0), 

only  the  fast  magnetoacoustic  mode  propagates.  Its  dispersion  relation  is 

2  2  % 

exactly  05  =  k+(s  +  a  )”*  and  v+  is  parallel  to  k+, 

v+  =  v+(0,  sin  0+,  cos  0+)  .  (4.' 

When  B  lies  in  the  plane  of  incidence,  v,  and  v  also  lie  in  this 
°  +  " 

plane,  while  v^  is  perpendicular  to  it.  Let  <p  and  a  be  the  angles  made  by  BQ 

***"  *  /VVW 

and  v,  respectively,  with  the  z  axis. 


Bo  =  sin  cp’  cos  4*) 


v«*  =  v0,(  1,  0,  0) 


v<*  - 


(4.9) 


where 


Vj.  =  v+(0,  sin(3  +  ,  cos(3+) 


(4.10) 


(?+  =  9  +  n/2  -  (s  /2a2)  sin|^2(0+  -  <p)J  +  0(s/a)4  ,  (4.1l) 


£3.  =  <p  -  ( s^/2a2 )  sin  Q(0_  -  cp)J  +  0(s/a)1 


(4.12) 


A  A  A  A 

Thus,  to  zeroth  order  in  s/a,  v+  is  perpendicular  to  BQ  and  v_  is  parallel  to  BQ. 

(MM 

The  energy  transport  vector  is  given  by  Eq.  (2.26)  which,  when  written 


here  in  terms  of  the  fluid  velocity,  v,  is 
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T  *(a2p( 


o/0))|  k(vxB0)2  -  (vxBq)  [k*(vxB0)|  +  (a2/a2)  v(k-v)| 

*•»  **s,  |  I  ****  **■»  «w  1 


(4.13) 


For  the  different  modes,  this  reduces  to 


To<  =  B0P0a(v«)2a"(k<»«'B0' 


and 


T±  "(a2p0/o))J  k^( v+xB0 

***-  i  • — -  ~~ 


)*  +  (  s2 / a2 ) 


1  v+(k+’v+M 

^  1 


(4.14) 


(4.15) 


whereo-(f)  =  f/|f|.  Equation  (4.14)  shows  that  the  energy  in  the  Alfven  mode 
propagates  in  the  direction  of  the  constant  magnetic  field. 

When  Bo  is  normal  to  the  plane  of  incidence,  the  relevant  field  quan¬ 


tities,  written  in  terms  of  v+,  are 


E+  =  (B0/c)  v+(°>  -cos  0+>  8ln  e+) 

(4.16) 

and 

B+  =  (B0k+v+/u5)  (1.  0,  0) 

(4.17) 

From  Eq.  (4.15) 

T+  -  (poa2^+/u))  (v+)2[^  +  (s2/a2)J 

(4.18) 

which  shows  that  the  energy  transport  is  entirely  in  the  k  direction. 

When  Bq  lies  in  the  plane  of  incidence,  the  field  quantities,  written 
in  terms  of  v,  are 


E  .  =  (B  /c)vw(0,  cos  Cp,  -sin  cp)  , 

(4.19) 

wJ?  o 

E+  =  -(Bq/c)v+  sin(p+  -  cp)  (1,  0,  0) 

(4.20) 

<V«M 

=  -(Bq/u))^^  008(6^  -  cp)  (1,  0,  0) 

(4.21) 

/vw* 
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and 


From  Eqs.  (4.14)  and  (4.15) 


-  Cp)  (0,  -cos  e+,  sin  e+) 

(4.22) 

cpyj 

(4.25) 

fk,  +  ?+  0(s2/a2f| 

(4.24) 

cos ( 0  -  9)  +  k  0(s2/a2)~l  . 

vwT  _ | 

(4.25) 

The  energy  in  the  fast  wave  of  the  magnetoacoustic  mode  propagates  mainly  in 

O  O 

the  k,  direction  with  a  component  0(s  /a  )  less  in  the  v+  direction,  while  the 

ivw  **** 

energy  in  the  slow  wave  of  the  magnetoacoustic  mode  propagates  mainly  in  the 

B  direction  with  a  component  0(s2/a2)  less  in  the  k_  direction. 

— 

5 .  PROPERTIES  OF  SOUND  AND  ELECTROMAGNETIC  WAVES 
Equations  (5.9)  and  (5.10)  give 


92v/3t2  =  s2r(v-v) 


(5.D 


A  A 

which  is  the  equation  of  a  longitudinal  wave,  i.e.,  vg  =  k8,  where  the  sub¬ 
script  s  denotes  sound  waves.  The  dispersion  relation  is  to  =  kss,  and  the 
energy  transport  vector  is 


Ts  «  pvs  =  (s2p0/u>)  ks(v8)2  (5.2) 

Equations  (5.11)  -  (5.14)  describe  electromagnetic  waves.  The  dis¬ 
persion  relation  is  ix>  =  kec,  where  the  subscript  c  denotes  electromagnetic 
waves.  With  ke  =  (0,  sin  8  ,  cos  @e),  the  electromagnetic  modes  may  be  de- 


-l6- 


scribed  completely  by  specifying  Ex  and  Ey.  Then  Ez 


-Eytan  6e>  Bx  -  -Eysec  0e, 


By  ■  Ex  cos  0e,  Bz  =  -Ex  sin  0 e>  and 


(c/M  ke[Ex 

■wx  I 


2  +  E  ^  sec2  9 


■] 


(5.3) 


6.  CALCULATIONS  FOR  BQ  NORMAL  TO  THE  PLANE  OF  INCIDENCE 


When  B0  is  normal  to  the  plane  of  incidence,  only  the  fast  magneto- 
***** 

acoustic  mode  propagates.  The  requirement  that  the  boundary  conditions,  Eqs. 
(2.17),  (2.18),  and  (2.21),  be  valid  for  all  times  over  the  entire  interface, 
in  the  usual  manner,  leads  to  the  conditions  that  a>  and *nxk  are  constant. 

*M>  MS 

Y 

These,  along  with  the  dispersion  relations,  are  sufficient  to  calculate  k 
and  kc  in  terms  of  .  One  finds,  then,  that  the  angle  of  incidence  equals 

MW 

the  angle  of  reflection. 


e+r 

=  Jt  -  0+1 

(6.1) 

V 

II 

+*■ 

P* 

(6.2) 

keC 

-  k^s2  +  a2)Vc 

(6.3) 

sin  9eL 

=  c(sc  +  a^)  ^  sin  0+*  , 

(6.4) 

C 

=  k+  (s  +  a  )  /s 

(6.5) 

sin  ©g1" 

2  2-%  i 

=  8(8  +  a  )  8  in 

(6.6) 

Equation  (6.4)  shows  that  whenever  0+*  is  greater  than  ( s^  +  a^)^/c,  sin  0^ 
is  greater  than  unity  so  that  cos  0eC  is  imaginary.  In  this  case,  the  electro¬ 
magnetic  wave  is  exponentially  damped  in  the  z  direction.  Thus,  the  incident 
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fast  magnetoacoustic  mode  must  have  its  wave  vector  lie  within  a  small  cone 

2  2  h 

of  half  angle  (s  +  a  P/c  with  cone  axis  along  the  normal  to  the  boundary 
in  order  for  the  electromagnetic  wave  to  propagate  unattenuated.  When  the 
incident  wave  vector  lies  outside  this  cone,  the  transmitted  electromagnetic 
wave  is  exponentially  damped,  proportional  to  exp(-z/d)  where,  from  Eq,(3.17).* 
d*  (s2  +  a2)V|  to  sin  0+*  I  •  Equation  (6.6)  shows  that  when  s/a  «  1,  the  fast 
magnetoacoustic  mode  incident  on  the  boundary  gives  rise  to  a  sound  wave  ly¬ 
ing  within  a  cone  of  half  angle  (s/a)  with  cone  axis  along  the  boundary  nor¬ 
mal.  That  is,  a  fast  magnetoacoustic  wave,  in  this  approximation,  gives  rise 
to  a  sound  wave  which  propagates  at  a  small  angle  to  the  normal  regardless  of 
the  angle  of  incidence. 

By  applying  the  boundary  conditions,  Eqs.  (2.17),  (2.18),  (2.21), 
and  (4.l6)  -  (4.18),  one  finds 

-(Bo/c)(v+i  cos  e+i  +  v+r  cos  ®+r)  “  Eyt  •  (6.7) 

v+i  cos  0+*  +  v+r  cos  0+r  ■»  v8t  cos  08t  (6.8) 

(s2po/o)(x+iv+i  +  k+rv+r )  +  (a2p0/u)Xk+1v+i  +  k+rv+r) 

=  (s2p0/u))(kstv8t)  -  (B0/43t)(Eyt/cos  0^)  (6.9) 

From  the  latter  equations,  with  the  aid  of  Eqs.  (6.1)  -  (6.6),  one  can  show  that 
Eyt  /E+;y  -  (V+*  -  v+r)/v+i 

=  2c(s^  +  a2)^|^c(s‘  +  a2)15  +  a2cos  0+i/cos  0et  +  sc  cos  0+i/cos  08t 

=  2  (s  «  a  «  c)  j  (6.10) 
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dC 


b>+;x  =  [(»: 


a  t 

<  n  •  T„C> 

—  -oy 
A  , 

<  n  •  T^_  A> 


/N  r 

<  n,Tat> 

*— ■  ~Jg 


v,uo  u  |  f  a  v-ua 

•2  +  a2)^/cjfc 


cos  0+i/cos  0e^j  |jv+i 


=  2a  cos  6+^/c  cos  ©g*-  (s 

4a2c(cos  0+i/cos  0et)(s2  +  a2)' 


«  a  «  c) 

,2^ 


(6.12) 


e+ _ 

©g*-  +  sc  cos  0+ 


M-ct  C7_|.  /CUB  Ug 

:(s2  +  a2)^  +  a2cos  0+^/cos  0e' 

( s  «  a  «  c ) 


0+*/cob  08t^J^ 


^/cos  ©g1") 


(4a/c)(cos  0+' 

o 

4sc  (cos  0+ 


(6.13) 


<  n.T,  *>  c(  s2  +  a‘ 


_i/ cos  9bL)(b2  +  a2)^ 

•2>*  +  a2cos  0+^/cos  0 e*-  +  sc  cos  0+Vcos  0Bt 


=  (4s/a)(cos  0+t/cos  08fc)  (s  «  a  «  c)  , 


(6. 


Equation  (6.13)  gives  the  energy  coupling  coefficient  for  electromagnetic  to 
magnetohydrodynamlc  waves,  and  Eq .  (6.14)  gives  the  energy  coupling  coeffi¬ 
cient  for  acoustic  to  magnetohydrodynamlc  waves.  In  the  limit  of  8  *=  0, 

Eqs .  (6.10)  -  (6.12)  reproduce  Eqs.  (  4l  )  -  ( 43  )  of  Ullah  and  Kahalas^. 

7.  CALCULATIONS  FOR'  B0  IN  THE  PLANE  OF  INCIDENCE 
In  contrast  to  the  case  of  ^Bg  perpendicular  to  the  plane  of  Incidence 
where  only  the  fast  magnetoacoustic  mode  can  propagate.  In  general  now  all 
three  magnetohydrodynamlc  modes  can  propagate.  It  will  be  shown  in  this  sec¬ 
tion  that  the  two  magnetoacoustic  modes  are  coupled,  but  propagate  Indepen¬ 
dently  from  the  Alfv/n  mode.  Before  discussing  the  coupling  of  wave  ampli¬ 
tudes  across  the  boundary,  we  shall  consider  the  relations  between  the  inci- 
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dent  wave  vector  and  the  reflected  and  transmitted  ones. 


A.  Relations  between  Wave  Vectors 

The  calculation  of  the  reflected  and  transmitted  wave  vectors  in  terms 
of  the  incident  wave  vector  involves  the  requirement  that  the  boundary  condi¬ 
tions,  Eqs.  (2.17).  (2.l8),  (2.23),  and  {2.2.h\  are  satisfied  over  the  entire 

A 

boundary  plane  for  all  times.  This  leads  to  the  statement  that  o>  and  nxk 
are  constant.  In  addition,  one  must  impose  the  requirement  that  an  incident 
wave  propagating  in  the  plasma  transports  energy  towards  the  boundary  and  a 
reflected  wave  propagating  in  the  plasma  transports  energy  away  from  the  bound¬ 
ary.  These  requirements  must  be  kept  in  mind  because  it  is  quite  possible  for 
/ 

the  Alfven  and  slow  magnetoacoustic  modes  to  transport  energy  in  such  a  manner 

A 

that  T-n,  the  projection  of  the  energy  transport  vector  in  the  z  direction,  has 

A 

opposite  sign  tojoru  This  situation  will  be  discussed  in  more  detail  below. 

For  B0  lying  in  the  plane  of  incidence,  there  are  several  different 
situations  that  can  arise.  These  are  specified  in  terms  of  the  wave  mode  in¬ 
cident  . 


l)  Alfven  Mode  Incident 

Consider  now  the  relation  between  wave  vectors  of  the  reflected  Alf- 
/  / 

ven  mode  and  the  incident  Alfven  mode.  Since  to  is  constant, 


Ka  IcosCe^1  -  qp)  |  =  k^  |  cos(0o(r  -  Cp)  |  .  (7.1) 

A 

Since  nxk  is  constant. 

k^,1  sin  =  kj0  sin  .  (7.2) 
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From  these  two  equations  one  finds  two  solutions 


i  r 

(cot  0^  +  tan  cp)  =  +  (cot  6^  +  tan  cp) 


(7.3) 


From  Eq.  (4.23)  <  n-To<>  =  pQa  <  (v^)^  >  jcos  q>  cos(0»f-  cp)|or-(cos  cp) o"jj:08( 0^ 

-  <P)]  .  In  order  for  the  Incident  wave  energy  to  impinge  on  the  boundary  plane, 
Cjr(cos  cp)o"||cos( 8*t  -  cp]  >  0,  and  for  the  reflected  wave  energy  to  be  propa¬ 


gated  into  the  plasma  away  from  the  boundary ,  CT"  ( cos  cp)a"  ^cos(ecAr  -  ,] 
Thus  it  is  required  that 


<  0. 


cr  ( cos  cp)  =  cr 


[cos(^  -  cp]  =  -cr|^os(e0,r  -  cp], 


(7.4) 


Equation  (7-4)  is  consistent  with  Eq.  (7.3)  only  when  the  minus  sign  is  chosen. 
This  leads  to 


cot  ej?  = 


cot  ^  -  2  tan  cp. 


(7.5) 


Figure  2  shows  a  plot  of  Fn.  (7.5)  for  various  values  of  qp  in  the  first  and 
fourth  quadrants.  0^  is  restricted  without  loss  of  generality  to  lie  in  the 
range  0  <  0*"  <  it.  Only  for  cp  =  0  or  for  very  small  angles  of  incidence  does 
the  angle  of  incidence  equal  the  angle  of  reflection  (i.e.,  01  =  n  -  0r). 

A  graphical  construction,  Fig,  3>  is  valuable  in  seeing  the  relation¬ 
ship  between  the  reflected  k«*r  and  the  incident  k^.  The  dot-dashed  line  is 

A  A  i 

the  locus  of  all  vectors  k  for  which  nxk  =  nxk^  .  The  dashed  lines  are  the 

A  A 

locus  of  all  vectors  k  for  which  | k • B  |  »  |  k^1  • B0 1  .  In  Fig.  3a.  the  inter- 
sections  of  the  dashed  lines  with  the  dot-dashed  line  gives  both  the  incident 
wave  vector  and  the  single  possible  reflected  wave  vector  at  A.  In  Fig.  34, 

BQ  lies  along  the  y  axis,  and  the  two  equations,  (7.1)  and  (7.2),  provide  the 
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same  information.  In  this  case,  there  is  no  unique  solution.  In  Fig.  3c> 

r  A 

both  intersections  lie  in  the  upper  half  plane.  Note  that  <  T™  ■ n  >  <  0,  as 

r  * 

required,  whereas  k^  -n  >  0. 

t  i  A 

The  relationship  between  k„  and  ko,  (obtained  from  cjd  and  nxk  con¬ 


stant)  is  given  by 


iin  =  (c/a)  ^sin  jcos(0o<i  -  cp)[j 


ket  =  kc*a|  cos(eJ‘  -  CP)|/C 


(7.7) 


For  the  electromagnetic  wave  to  propagate  unattenuated,  sin  0ec  must  be  less 
than  or  equal  to  unity.  Equation  (7.6)  shows  that  the  incident  wave  vector 
must  lie  within  a  cone  (the  cone  for  radiation)  whose  axis  lies  along  the  nor¬ 
mal  to  the  boundary  and  with  half  angle  for  a/c  «  1  given  by  ajcos  cp|/c. 

When  the  electromagnetic  wave  is  damped,  proportional  to  exp(-z/d),  d  is  given 
by  Eq.  (3.17) 


d  —  a  |  cos(0e/  -  cp)  |  /  a> | sin  ^  j 
2)  Fast  Mode  Incident 


(7.8) 


Consider  first  the  relation  between  wave  vectors  of  the  reflected  and 
the  incident  fast  magnetoacoustic  mode.  To  lowest  order  in  s/a,  the  dispersion 
relation  of  the  fast  magnetoacoustic  mode  is  independent  of  the  orientation 
of  k+.  Since  co  is  constant 


k  1  =  k  r 

K+  K+ 


and  since  nxk  is  constant, 


k+^  sin  0_|/  =  k+r  sin  0+r 


(7.10) 
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Thus,  the  angle  of  Incidence  equals  the  angle  of  reflection, 


.  r  i 

e+  =  *  -  e+ 


(7.11) 


From  the  requirement  that  <  T  On  >  >0  and  <  T+r.jn__>  <  0  and  from  the  fact 
that  Tj_  is  in  the  direction  of  k  to  lowest  order  in  s/a,  0  ^  0+  ^  Jt/2  ^  0 


Next  the  relation  between  wave  vectors  of  the  reflected  slow  magneto¬ 
acoustic  mode  to  the  incident  fast  magnetoacoustic  mode  is  considered.  With 
03  constant 


k^  a 


k_r  s  | cob(0  r  -  cp)  | 


(7.12) 


and  nxk  constant 


i  i  r  r 

k+  sin  6+  =  k_  sin  0_ 


(7.13) 


From  these  two  equations  one  finds  two  solutions 

a/ |^s  sin  0+^  jcos  cp  j^J  =  +  (cot  0  r  -4-  tan  cp)  .  (7.1*0 

Without  loss  of  generality  0  is  taken  so  that  0  8  ^  jt  for  all  modes.  From 

Eq.  (4.25)  <  ri*T_  >  =  (p0s^k_/u>)<  (v_)^  >  |cos  qp  cos(0_  -  qp)  J  cr*  (cos  cp)o_|cos(0_ 

-  cp)  I  .  The  requirement  that  <  n-T  >  <0  leads  to 

_4  -w 

0~  (cos  cp)  =  -<rJj;os  (6-r  *  <?)]  .  (7.15) 

Since  0~(  s  in  0_r)  *>  +  1,  there  results  a~  jj:os  (S_r  -  cp)]  -  a~  (cos  9)0"  (cot  0_r 
-f-  tan  cp) .  This  combined  with  Eq.  (7.15)  leads  to  the  condition  that  O”  ( cot  0_r 
•+.  tan  cp)  =  -  1.  Thus  Eqs.  (7.14)  and  (7. 15)  are  consistent  only  when  the  minus 


-23- 


sign  is  chosen,  which  gives 


cot  9  r  =  -Qa/(s  sin  0+*  |cos  q)|)J  -  tan  cp 


(7.16) 


Equation  (7.16)  shows  that  for  any  angle  of  incidence  the  reflected  wave  lies 
within  a  cone  of  half  angle  sjcos  cpj/a  whose  axis  lies  along  the  normal  to  the 
boundary . 

Figure  4  shows  the  construction  by  which  k_r  is  found  given  k^ , 

A  A  . 

Again,  the  dot-dashed  line  is  the  locus  of  all  vectors  k  for  which  nxk  =  nxk,  ■L. 

The  circle  has  a  radius  ak  Vs,  and  the  dashed  lines  are  the  locus  of  all  vec- 
.  A  t 

tors  k  for  which  k- B  =  ak  /s.  In  Fig.  4a,  the  intersection  of  these  loci 
I  —  ^2  1  + 

occurs  at  A  and  A' ,  By  considering  all  possible  orientations  of  B0,  one  can 

show  that  k_r  is  given  by  the  intersection  at  A  in  the  lower  half  plane  and 

that  the  k_r  always  lies  in  the  lower  half  plane.  Figure  4b  shows  the  case 

when  cp  -  +  n/2.  In  this  case,  there  is  no  solution. 

The  relations  between  the  transmitted  wave  vectors  and  the  wave  vec- 

A 

tor  of  the  incident  magnetoacoustic  mode  (given  by  a)  and  nxk  constant)  are 

=  k+1  aj z  (7.17) 


and 


r  i  1 

sin  6  =  c  sin  6  / a 

e  +  l 


for  the  electromagnetic  wave  and 


s'  =  k+'  a/8 


(7.18) 


(7.19) 


and 


sin  0„ 


s  sin  & 


Va 


s  + 

for  the  sound  wave.  For  the  electromagnetic  wave  to  propagate  unattenuated 


(7.20) 
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sin  0  must  be  less  than  or  equal  to  unity.  Equstion  (7.18)  shows  thst  the 
incident  wave  vector  must  lie  within  a  cone  whose  axis  lies  along  the  normal 
to  the  boundary  and  with  half  angle  sin  ^(a/c).  The  angle  at  which  the  sound 
wave  is  transmitted  is,  from  Eq.  (7.20),  close  to  the  normal  for  any  0+ *  when 
s/a  «  1.  When  the  electromagnetic  wave  is  exponentially  dstnped,  proportional 
to  exp(-z/d),  d  is  given  by  Eq.  (3.17) 

d  —  a/  cujain  j  (7.21) 

3)  Slow  Mode  Incident 

Consider  first  the  relation  between  wave  vectors  of  the  reflected 
fast  magnetoacoustic  mode  to  the  incident  slow  magnetoscoustic  mode.  Since 
O)  is  constant, 

k+ra  =  k.^s  jcos(6.i  -  <p)|  ,  (7.22) 

/\ 

and  since  nxk  is  constant, 

k+r  sin  6+V  =  k.1  sin  e.1  .  (7.23) 


These  two  equations  have  two  solutions 


a/(s  sin  6+v  |cos  ®|  )  “  +  (cot  +  tan  (p) 


(7.24) 


Without  loss  of  generality  6  is  taken  so  thst  0  0_  <n.  From  Eq.  (4.25) 

<Jn-T_  >  -(pos2k./oj)<  (v_)2  >  |  cos  cp  cos(0_  -  (p)|  or  (cos  cp)o”Jjcos  (6.  *  9)]  • 

The  requirement  that  >  0  leads  to 

a~  (cos  q>)  ■  o"  ^cos  (e.1  -  <p)J  (7.25) 
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With  0_^  in  the  first  and  second  quadrants,  cr  ^cos  (e.1  -  <p)^j  -  cr(cos  cp) 
a( cot  0_*  +  tan  cp) .  The  latter  equation  combined  with  Eq.  (7.25)  leads  to 
the  condition  that  <j(cot  0.^  +  tan  cp)  =  +  1.  Thus  Eqs.  (7.24)  and  (7.25) 
are  consistent  only  when  the  plus  sign  is  chosen  which  gives 


sin  0+r  »  (a/s)  £  |  cos  9  j  (cot  0_^  +  tan  <p)J 


(7.26) 


For  the  reflected  fast  magnetoacoustic  mode  to  propagate  unattenuated,  cot  0_ 
>,  (a/s  |  cos  cpj  )  -  tan  (p,  which  may  be  approximated  by  0_^  ^  (s/a)  |  cos  9  |  . 
That  is,  the  incident  slow  magnetoacoustic  mode  must  be  within  a  cone  of  half 
angle  (s/a)  |  cos  9  |  .  When  the  fast  magnetoacoustic  mode  is  attenuated,  pro¬ 
portional  to  exp(-|z|/d),  d  is  given  by  Eq.  (3.17) 


da^-  s  |cos(0_i  -  9)  |  j  o>|sin  0 _^| 


(7.27) 


Figure  5  shows  the  construction  by  which  is  found,  given  k_  . 

A  A  t 

Again,  the  dot-dashed  line  is  the  locus  of  all  vectors  k  for  which  nxk  “  nxk 

A 

The  circle  has  a  radius  (s/a)  |k_i-Bo|»  ca/a.  Then  k+r  is  determined  by  the 

- -  MM 

intersection  of  the  circle  and  the  dot-dashed  line.  There  can  be  two  inter¬ 
sections,  and  the  one  in  the  lower  half  plane  is  chosen,  according  to  the  re¬ 
quirement  ^  T.r  -n  ^  <  0,  since  the  fast  magnetoacoustic  mode  transports  en- 

Hk  **** 

ergy  in  the  direction  of  k+r.  In  Fig.  5a  the  intersections  occur  at  A  and 
A' ,  and  the  one  in  the  lower  half  plane  is  chosen.  Figure  5b  shows  a  con¬ 
struction  in  which  no  solution  exists.  In  this  case  sin  0+r  >  1.,  and  Eq. 
(7.26)  leads  to  cot  0_^  <  (a/s  |  cos  9  j)  -  tan  9.  Figure  5C  gives  the  con¬ 
struction  for  9  “  +  «/2  for  s/a  <  1.  In  thi3  case  no  solution  exists. 
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Next,  the  relation  between  wave  vectors  of  the  reflected  slow  mag¬ 
netoacoustic  mode  to  the  incident  slow  magnetoacoustic  mode  is  considered. 
With  co  constant. 


|  cos(e_1  -  cp)  |  =  k_r  |  cos(6_r  -  cp)|  , 


(7.28) 


and  with  nxk  constant, 


k  1  sin  fi  1  =  k  r  sin  0  r 


(7.29) 


These  two  equations  are  the  same  as  the  ones  for  Alfven  mode  incident  and  re¬ 
flected.  The  discussion  is  the  same  as  for  that  case,  with  the  resulting 
equation 


cot  0_r  *>  -  cot  S,1  -  2  tan  cp 


(7.30) 


Figure  2  shows  the  relation  between  the  angles  of  incidence  and  reflection, 

r  i 

and  Fig.  3  allows  the  construction  of  k_  from  k_  . 

The  relations  between  the  transmitted  wave  vectors  and  the  wave  vec¬ 
tor  of  the  incident  slow  magnetoacoustic  mode  (given  by  CD  and  nxk  constant) 


k  -  k.1  |cos(0.i  -  cp)  I  s/c 


(7.31) 


dj"  -(c/s^sin  e^j  |  cos(e.i  -  cp) 


(7.32) 


for  the  electromagnetic  wave,  and 


k8t  “  k.1  |  cos(e.1  -  CP)  | 


(7.33) 
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and 


(7-34) 


sin  0gi  »  sin  ©_*■/  |cos(0_^  -  cp)| 

for  the  sound  wave.  Equation  (7-34)  is  shown  graphically  in  Fig.  6.  Equation 
(7.32)  shows  that  for  the  electromagnetic  wave  to  propagate  unattenuated,  the 
incident  wave  vector  of  the  slow  magnetoacoustic  mode  must  lie  within  a  cone 
whose  axis  lies  along  the  normal  to  the  boundary  and  with  half  angle  given  by 
s  j  cos  <p  |  /c.  When  the  electromagnetic  wave  is  exponentially  damped,  propor¬ 
tional  to  exp(-z/d),  d  is  given  by  Eq.  (3. 17) 


d—  s  jcos(0_i  -  cp)  [/  cn|sin  ©_ .  (7.35) 

Equation  (7-34)  shows  that  the  sound  wave  is  damped  unless 


cot 


6.1  ^,(1/  |  cos  cp  |  )  -  tan  cp 


(7.36) 


Whenever  -Jt/2-$  <p  ^  Jt/2,  Eq.  (7-36)  is  equivalent  to  0_^  ^  cp/2  +  n/4.  When¬ 
ever  «/2  ^  <P-^3fl/2,  Eq.  (7.36)  is  equivalent  to  0_^  cp/2  -  «/4.  When  the 
sound  wave  is  exponentially  damped,  proportional  to  exp(-z/d),  d  is  given  by 

Eq.  (3.17) 


d — 


s  jcos( 0.  ^ 


cos  q>  cos(2 0_ 


i 


(7.37) 


B.  Relations  between  Wave  Amplitudes 
The  boundary  conditions,  Eqs .  (2.17),  (2.18),  (2.23),  and  (2.24),  are 
now  applied,  along  with  the  equations  of  Sections  4  and  5  in  which  the  field 
vectors  are  expressed  in  terms  of  specific  amplitudes  and  angles.  The  re¬ 
sulting  equations  separate  into  two  sets.  The  first  set  contains  only  v^, 
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v  r,  and  E  fc:  the  second  contains  only  v  *  v  r,  E  and  v  Thus,  an  in- 
cident  Alfven  mode  gives  rise  only  to  a  reflected  Alfven  mode  and  a  trans¬ 
mitted  electromagnetic  wave.  An  incident  magnetoacoustic  mode  gives  rise  only 
to  reflected  magnetoacoustic  modes  and  transmitted  electromagnetic  and  sound 
modes . 

.  ✓ 

1)  Alfven  Mode  Incident 


The  equations  relating  amplitudes  on  the  two  sides  of  the  boundary 


BQ  cos  cp(vQ1  +  var)/c  =  Eyfc 


(7.38) 


(B_/u)|  k  *v  cos(6  ^  -  cp)  +  k  rv  r  cos (0  r  -  9)  I  =  E  t/cos  0  t  (7*39) 
v  0  a  a  x  a  a  a  a  I  y  e  v  ' 

These  equations  determine  var  and  EyC  in  terms  of  vai .  From  Eqs.  (7-38)  and 
(7-39).  the  ratios  of  field  quantities  are 


E  /E  ■  (v  ^  +  v  r)/v  t 

y  /c,y  v  a  a  >'  a 


2  cos  0, 


•Vtcos  6et  +  (a/c)  j  COS  <pfj 


(a/c  «  1)  ,  (7.40) 


Ezt/Ez  1  -  (tan  0et/tan  q>)(Eyt/Eyi) 


2  sin  9^  cot  <p/  £cos  0g  +  (a/c)  |  cos 


■  2  tan  6  cot  <p 


(a/c  «  1)  ,  (7.41) 


B^/B^  “  (2a/c)  I  cos  q>  |  /[cos  6et  +  (a/c)  lcos  *f] 


(2a/c)  |  cos  q>  |  / cos  0 


(a/c  «  1)  .  (7.42) 
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The  expressions  for  the  coupling  coefficients  are 


<  n-T  c  >  (4a/c)|cos  9  |  cos  6  1 


^  ,  Q 

<  n'T.,1  >  (cos  8  c  +  a|  cos  9l/c) 

=  (4a/c)|cos  9l'/cos6et  (a  «  c)  ,  (7*45) 


<  n-T  r  >  fcos  &J1  -  a|cos  9  |/< 


(a  «  c)  .  (7-44) 


<  n-T^i  ^  cos  Sg1"  +  a|cos  9l/c 


2)  Fast  or  Slow  Mode  Incident 

The  equations  relating  the  amplitudes  on  the  two  sides  of  the  bound¬ 


ary  are 


(B0/c)^+1sin(ft+1  -  9)  +  v_isin(p_i  -  9)  +  v+rsin(p+r  -  9) 

+  v.rsin(p.r  -  9)J  =  ~ 


(7.45) 


(l^j /o)|k+^v+1‘sin(p+1‘  -  9)cos  6+^  +  k_  ^"v_  isin(  -  9)cos  6_^  +  k+rv+r- 

•  sin(p+r  -  9)cos  6+r  +  k.rv.rsin(p.r  -  9)cos  6_r^  “  '  cos  *  (7.46) 

k+1v+icos(  -  p+*-)  +  k_*-v_*-cos(6_i  -  p_  *■)  +  k+rv+rcos(  6+r  -  |3+r) 

+  k.rv.rcos(0.r  -  p_r)  -  kgtv8c  ,  (7.47) 


v+*-cos  +  v_^cos  +  v+rcos  p+r  +  v_rcos  B_r  •  vgtcos  8 gt  .  (7.48) 
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For  the  fast  magnetoacoustic  mode  incident  v_i  **  0  and  cn  is  given  by  Eq.  (4.5); 
for  the  slow  magnetoacoustic  mode  incident  =  0  and  oe  is  given  by  Eq.  (4.6). 

The  solution  of  Eqs.  (7*45)  •  (7*48)  for  the  fast  wave  incident  gives 
for  the  ratio  of  field  amplitudes 


and 


E^/E*1  ■=  (v*1  +  V^/v/ 

=  2  cos  e+  6+1  +  cos  ®efcJ 

(a  «  c)  , 


=  2 


Bz'/Bz1  = 


2  cos  JTos  0+*  +  (a/c)  cos  0^^ 


=  2 


B  t/B  1  -  (v  1  -  v  r)/v  1 
y  y  +  +  '  + 


(7.49) 


(a  «  c)  ,  (7.50) 


=  (2a/c)  cos  eey  |Tos  0+*  +  (a/c)  cos  ©et~j 
=  (2a/c)(cos  0et/cos  S^.1)  (a  «  c)  4  (7-51) 

The  expressions  for  the  coupling  coefficients  are 

<  n-T  t  >  (4a/c)  cos  6,  cos  0  t 

<  n-T  >  f"cos  9  1  +  (a/c)  cos  6  t"| 

~  L  +  e  J 

-  (4a/.c)(  cos  0et/cos  e^.1)  (a  «  c)  ,  (7-52) 
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A  t 

<  n-V  > 


(4s/a)  sin^cp  cos  cos 


<  r«T+^  >  cos  6+i  +  (a/c)  cos  ®e^|  j  cos  <p  |  +  cos  6at~j 


(7-53) 


In  the  above  only  the  lowest  order  In  s/a  has  been  kept.  Equation  (7*53)  gives 
the  coupling  coefficient  of  sound  waves  to  magnetohydrodynamic  waves.  It  is 
of  the  order  of  s/a  which  has  been  assumed  to  be  small  in  this  paper.  In  the 
lower  ionosphere  where  collisions  between  heavy  ions  and  neutral  particles 
become  important  (lower  E  region),  the  phase  velocity  of  magnetohydrodynamic 
wave  propagation  can  be  reduced  to  a  small  fraction  of  the  Alfven  velocity 
and,  in  fact,  can  become  comparable  with  the  sound  velocity.  In  this  case 
there  may  be  a  large  coupling  between  magnetohydrodynamic  and  acoustic  waves, 
although  the  present  model  is  not  adequate  to  treat  that  situation  because 
only  a  fully  ionized  plasma  has  been  discussed  here. 

For  s/a  «  1,  sin  Sgt  is  small  and  cos  8^—  1.  For  sin  S+4  < 
g/c,  the  electromagnetic  wave  is  undamped  and  Eq.  (7.53)  can  be  written 


A 

<  n-T  c  > 

Ji 

1  > 

+ 


(4s/a)  tan2(cp/2) 


=  (4s/a)  cot"(cp/2) 


(  -«/2  ^  cp  ^  n/2) 


(  «/2  ^  cp  ^  3*/2) 


(7-54) 


The  solution  of  Eqs.  (7-45)  -  (7-48)  for  the  slow  wave  incident  gives 
for  the  ratio  of  field  amplitudes 

Ext/Ex1  “  *(a/s)[^os  &.i  &“/|cos (S.1  -  cp ) Q  (7-55) 
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t 


I 


where 


r  r  ij 

sin(0_r  -  p)  cos  0_r  Icos  p|  -  cos  0gt  f  -1 

H  Jl  + - J - I -  [cos  0+r  -  (a/c)cos  0ecJ 

sln(0.^  -  p)  cos  0_^  |cos  pj  +  cos  Ogt 

V.  L  -y 

When  the  transmitted  electromagnetic  wave  is  undamped,  6_^-~  0,  Og*- 


—  0,  0.ri!  it,  and  0+r—  n.  Then  Eq.  (7*55)  becomes 


ExC/Exi-  (2a/s)/(l  +  |  cos  pj) 


(7.56) 


Byc/By1  =  -(a/c)  cos  0efc  6" 


(7-57) 


and  when  the  transmitted  electromagnetic  wave  is  undamped. 


Byt/By1  =  (2a/c)cos  0et|cos  p|/(l  +  | cos  p  |) 


(7.58) 


B^/B^  *  -(a/c)  sin  0gb  cot  0  *  6~ 


(7-59) 


which  becomes, when  the  electromagnetic  wave  is  undamped 


Bzc/Bzi  =  ( 2a/c )  sin  0ec  cot  0_1|cos  p|/(l  +  |cos  p  |)  }  (7.60) 

which  is  of  the  order  of  a/s  because  cot  0_^  is  of  the  order  of  c/s  when  0_* 
lies  within  the  cone  for  radiation  in  this  case. 


The  expressions  for  the  coupling  coefficients  are 


<(n_'Tet  >  (4s/c)  cos  0ec  sin2p  |cos  p| 

—  ~ "  r  ,  TP 

<  n.T-  >  11  +  j cos  plj 


(7.61) 
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and 


(7-62) 


Figure  7  is  a  plot  of 
magnetoacoustic  mode 
weakly  (of  the  order 


A 

<  «-Tst  > 


4|cos  cpj  cos  0gc 


<  "n.T_^  >  ^cos  cp  |  +  cos 


Eq.  (7.62).  As  one  might  expect}the  slow  wave  of  the 
couples  most  strongly  to  the  sound  wave  and  much  more 
of  s/c)  to  the  electromagnetic  wave. 
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APPENDIX 


In  References  l6  and  17,  "Coupling  of  Magne tohydrodynamic  to  Electro¬ 
magnetic  Waves  at  a  Plasma  Discontinuity.  I  and  II",  (hereafter  called  I  and 
II),  there  is  an  error  made  which  seriously  affects  some  of  the  results  of 
those  papers.  This  appendix  is  devoted  to  a  discussion  of  the  error  and  its 
correction.  The  work  of  Section  J  is  directly  related  to  this  discussion, 
since  the  problem  is  the  same  when  s  =  0. 

The  error  occurs  for  the  case  where  BQ  lies  in  the  plane  of  inci- 
/ 

dence  and  the  Alfven  mode  is  incident.  (In  I  and  II,  this  is  called  Case  lb, 

and  the  Alfven  mode  is  called  the  extraordinary  mode.)  The  error  is  that  the 

angle  of  reflection  was  assumed  to  be  equal  to  the  angle  of  incidence.  This 

assumption  is  correct  for  the  fast  magnetoacoustic  mode  incident  (called  the 

ordinary  wave  in  I  and  II),  so  that  all  the  results  of  Cases  la  and  2a  which 

/ 

refer  to  the  ordinary  mode  incident  are  correct.  For  the  Alfven  mode,  the 
relation  between  the  angle  of  incidence  and  reflection  is  given  in  Eq .  ( :f • 5 ) 
and  is  shown  graphically  in  Fig.  2.  It  is  seen  that  only  for  cp  -  0  or  for 
9  small  does  the  angle  of  incidence  equal  the  angle  of  reflection. 

Before  proceeding  to  the  discussion  of  the  revision  of  the  equa¬ 
tions,  we  wish  to  point  out  that,  in  spite  of  the  error  made,  the  results  of 
I  are  correct  even  for  Case  lb.  This  may  be  seen  from  the  following  reason¬ 
ing:  In  I  discussion  was  restricted  to  the  case  when  electromagnetic  waves 

in  the  vacuum  region  propagated  unattenuated.  This  restriction  implies  that 
the  wave  vector  of  the  incident  Alfven  mode  lies  within  the  cone  for  radia¬ 
tion,  i.e.,  that  sin  9  4  Vg  cos  <p/c  (Eq.  (3l)  of  I  in  the  notation  of  that 
paper).  When  the  angle  of  incidence  is  very  small,  it  is  seen  from  Eq.  (7.5) 
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that  for  all  cp,  0Q  -  n-  0Q  ;  i.e.,  the  angle  of  incidence  equals  the  angle 

of  reflection  if  the  angle  of  incidence  is  very  small.  Thus  the  assumption 
made  is  correct,  snd  sll  the  results  of  I  are  correct. 

Instead  of  Eq.  (25)  of  I  (in  the  notation  of  that  paper),  the  cor¬ 
rect  equation  is 


Ey1"  k  cos  S'  sin  9(yr  -  ?*■)  +  cos  0(1  -  tan  0  cot  ©r)| 
E  *  k  cos  6'(yr  sin  0  -  sin  0  cot  0f)  +  k'j 


(A.l) 


where  0r  is  the  angle  of  reflection  (as  defined  in  Fig.  1).  The  expression 
for  (Eq.  (20)  of  II)  is  given  correctly  in  that  paper  by 

i  _  cos ^(0  -  (p)sin^cp  +  x^[Tos^0  -  cos ^(0  -  cp)  -  sin^cp]  +  x^sin^0 

-cos ^(0  -  <p)sin  cp  cos  cp  +  x^(sin  q>  cos  cp  +  sin  0  cos  0)  -  x^sin  0  cos  0 


However,  the  correct  expression  for  y  is 


r  cos2(0  -  cp)sin2<p  +  x^Jjsin^©  cot ^9V  -  cos^(0  -  q))  -  sin^cpj  +  x^(l  -  sin^e  cot^0r) 

r  “  - p - r - 

-cos 2(0  -  q))sin  qi  cos  q)  +  x^lsin^©  cot  0r  +  sin  cp  cos  qn  -  x^sin^©  cot  ©r 


(A. 3) 


In  addition  Eq.  (7-5)  is  written  in  slightly  different  notation 


cot  ©r  m  -cot  0-2  tan  cp 


(A. 4) 


These  equations  are  adequate  for  discussing  the  behavior  of  E  First,  if 
there  were  a  transmission  resonance  in  Csse  lb  (i.e.,  Alfv^n  mode  incident) 
the  bracketed  psrt  of  the  first  term  in  the  denominator  of  Eq.  (A.l)  must 
vanish.  This  occurs  when  yr«  cot  ©r.  However,  yv  is  approximated  by  yr  5 
-  tsn  cp  which  lesds  to  the  solution,  using  Eq.  (A. 4),  ©  -  cp  •  +  «/2.  In  this 
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case,  the  Incident  Alfven  mode  does  not  propagate  sincejc  is  perpendicular 
to  Bq .  Thus  no  transmission  resonance  is  possible  in  Case  lb,  contrary  to 
the  results  of  II. 

Equation  (A.l)  will  now  be  evaluated.  The  left  hand  term  of  the  de 
nominator  of  Eq.  (A.l)  vanishes  only  when  the  angle  of  incidence  lies  on  the 
cone  for  radiation  (cos  0'  =>  0)  or  when  0  -  cp  =  +  n/2,  as  discussed  above. 
Excepting  these  cases,  the  right  hand  term  in  the  denominator  is  0(a/c)  with 
respect  to  the  left  hand  term  and  so  is  ignored.  Thus  Eq.  (A.l)  becomes 


sin  9(yr  -  y^)  +  cos  0(1  -  tan  0  cot  0r) 
(yr  -  cot  0r)sin  0 


(A. 5) 


On  using 


i 


y 


yr  =  -  tan  <p  together  with  Eq .  (A. 4),  one  gets  that 


E  t/E  i 

y  y 


(A.6) 


which  is  the  correct  answer  for  all  Cp.  The  only  remaining  point  is  to  see 
when  y  and  yr  may  be  approximated  by  -  tan  cp.  From  Eq.  (A. 2)  and  (A. 3),  one 
sees  that  the  approximation  is  permissible  when  the  terms  independent  of  x 
do  not  vanish.  When  these  terms  do  vanish  cp  “  0,  +  n/2,  and  9  -  cp  -  +  n/2. 

The  latter  condition  corresponds  to  the  case  of  no  propagation,  as  mentioned 
above;  cp  »  +  n/2  corresponds  to  an  indeterminate  situation  because  the  con- 

/N. 

ditions  cc  constant  and  nxk  constant  are  the  same  to  lowest  order  in  x.  If 
one  allows  finite  x,  then  the  angle  of  incidence  again  equals  the  angle  of  re¬ 
flection,  and  the  discussion  of  II  is  valid  (for  CP  ”  +  n/2).  For  CP  ”  0  also, 
the  discussion  of  II  is  correct. 
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Equation  (21)  of  II  is  not  changed,  since  it  has  no  dependence  on 


reflected  quantities.  Thus,  for  all  values  of  Cp,  (in  the  notation  of  II) 
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ANGLE  OF  REFLECTION  VS  ANGLE  OF  INCIDENCE  FOR  VARIOUS  ORIENTATIONS  OF  THE 

/ 

EXTERNAL  MAGNETIC  FIELD,  ALFVEN  MODE 
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■Fig.  6.  RELATION  BETWEEN  THE  ANGLE  OF  TRANSMISSION  OF  SOUND  WAVES  AND  THE  ANGLE 
OF  INCIDENCE  OF  THE  SLOW  WAVE  FOR  VARIOUS  ORIENTATIONS  OF  THE  CONSTANT 
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